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Semi-Annual Status Report - NSG 5014 - December 31, 1978 

This report covers activities of the Landsat Sensing 
Research Group (Earth Resources) and of the Volcanic Gas Sensing 
Research Group (Planetary Science) which work in collaboration 
with the Goddard Institute for Space Studies, New York, 

Dr. Robert Jastrow, Director and Technical Officer for this 
project. 

A. Landsat 

The Dartmouth Landsat Research Group continued application 
studies for Landsat data under the general category of analysis 
of vegetation cover, especially forestry and geobotany, that 
is, the effects of soil/earth mineral content on vegetation. 

1. Agricultural Crop Studies . No research activities in 
period. (No funds allocated.) 

2. Forestry (Emily Bryant). This work includes the effort 
of Ms. Bryant and her student assistants, in collaboration with 
A.G. Dodge, Area Forester, Cooperative Extension Service, 
University of New Hampshire, whose services are furnished under 
a task order to UNH by Dartmouth. 

We have continued in both lines of research mentioned 
in the previous report: clearcuts, consistency. 

A. Clearcuts — Field checked and took photos of 
recent clearcuts in northern New Hampshire in August. Photos 
used in a prototype "photo-key" — a visual aid to explain 
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computer classification of clearcuts to field personnel and 
researchers. Brought the clearcut classification to a point 
where we are ready to work with subsequent Landsat passes. 

B. Consistency -- Our assistant worked on getting 
consistent classification of an area in northern New Hampshire 
in three summer Landsat passes, 1973, 1975, and 1976. Pre- 
liminary results with the 1973 and 1975 data show acreages of 
major forest types (softwood, mixed wood, and hardwood) were 
within about 15 - '20% of each other on test areas of sizes 
ranging frcpm 20,000 to 400,000 acres. Checked our softwood 
training area to see if it is still valid to use with 1978 
data. 

Completed printout of Franklin County, Massachusetts, 

It showed more softwood than was indicated in a classification 
of the same area done by Bendix. 

We have also worked on a project to test the atmospheric 
simulations which Richard Kiang has developed. This involves 
classifying a 720,000 pixel (826,000 acre) area, and 
reclassifying it 6 times with 6 different atmospheric 
perturbation simulations. 

Spreading the word, maintaining contacts with remote 
sensing and forestry world. 

A. Spreading the Word — Dodge and/or Bryant gave talks/ 
lectures as follows: 

Squam Lake Science Center, Holderness, N.H. 
corporation meeting, July 22. 
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Meeting of the New England - St. Lawrence Valley 
Association of Geographers, October 20. 

Annual Meeting of all Cooperative Forest Management 
Supervisors in the Northeast (20 states) , 

September 19. 

Remote Sensing class, Dartmouth, November 14. 
Submitted two presentation summaries for poster 
session at 13th International Remote Sensing 
Symposium at ERIM, University of Michigan, April, 
1979 (Enclosures 1 & 2) . 

B. Contacts ' 

Went to Ottawa, Canada and visited Canada Centre 
for Remote Sensing and the Canadian Forest Fire 
Research Institute, September 5-7. Compared their 
techniques with ours using our computer compatible 
Landsat tapes. 

Went to meeting of the Remote Sensing Group of 
northern New England, September 28-29. Accepted 
some responsibility for the next annual program. 
Talked with Jack Estes of NASA Headquarters, 

October 1. 

Went to GISS and talked with people up from NASA; 
also exchanged ideas with people at GISS, 

November 8-9. 

Discussed Landsat data classification techniques 
with Ben Merembec from ORSER, Penn State. 

Maintained contact with Darrel Williams of GSFC 
by telephone. 
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3. Land Use . No research activities in period. (No 
funds allocated) . 

4. Water Quality . No research activities in period. (No 
funds allocated.) 

5. Mineral Resources . Work under this heading is 
undertaken by Professors Birnie and Stoiber, with Joseph 
Francica and student assistants. 

a) Airborne multispectral scanner experiments were 
undertaken in the fall of 1978 over the American River copper 
prospect in Washington State. The deposit is overlain by a 
dense fir forest with some tamarack. The purpose of these 
experiments is to see if the presence of the known mineral 
deposit is manifested asa geobotanically anomaly that can be 
detected by remotely obtained high resolution reflected radiance 
data. These data have been collected and are being analyzed 
now. 

The spectral characteristics of the fir forest 
will be compared with those of a lodgepole pine forest collected 
over the Heddleston porphyry copper molybendum deposit in 
Montana in 1975. A geobotanical anomaly was detected in the 
Heddleston data. If prospecting for ore deposits using airborne 
multispectral scanner data to detect geobotanical anomalies is 
to be universally applicable, the method must be tested at 
different types mineral deposits with different types of forest 


cover. 
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A classification algorithm for aircraft data has 
been developed. This classification discriminant is based on 
the observation in the Heddleston data that the mineralized 
spectra have higher reflected radiance values in the 0.67 pm 
region of the spectrum. 

b) Landsat Data. A study of the northwestern 
Adirondack Mountain region has been undertaken using digital 
processing of Landsat data. The Landsat data is being used to 
attempt to map a complex metamorphic area of marbles and 
gneisses. These rock types are manifested on the surface by 
different surface covers and vegetation patterns. A map has 
been procduced but it is being improved by more careful 
attention to classification parameters. The GISS "fixed 
volume" classification is being used. 

c) Application of Remote Sensing to Reconnaissance 
Geologic Mapping and Mineral Exploration. This paper by 
Prof. Birnie and Dr. Jon Dykstra (former graduate student) has 
been published in the Proceedings of the 12th International 
Symposium, Remote Sensing of the Environment (Enclosure 3, 
reprint) . 

B. Remote Sensing of Volcanic Emissions. This work 
was undertaken by Prof. Stoiber, Lawrence Malinconico, Technical 
Assistant Geoffrey Bratton and student assistants. 

1) Data from aircraft SC> 2 monitoring (COSPEC) 
in February 1978 was presented at a meeting of the GSA and of 
Canada (1978 Joint Meeting) in Toronto, October 28 and of the 
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American Geophysical Union in San Francisco in December 1978. 
This data is now being prepared for publication. (Abstracts 
at Enclosures 4 & 5) . 

2) Arrangements have been made to run another set 
of experiments on sulfur dioxide determination from the air 
similar to those of February 1978. This will be as before in 
cooperation with NCAR using the Queenaire in February 1980. 

3) The Remotely Piloted Aircraft project was 
inactive because plans to visit, for a field trial, the optimum 
locality at Masaya, Nicaragua, were cancelled due to extreme 
political unrest in Nicaragua. Samples were collected on the 
ground at the Volcano San Miguel in El Salvador. These were 
used to further refine our analytical techniques for the 
determination of chlorine and sulfur dioxide in gases collected 
either by hand or by remotely piloted vehicles at fuming 
volcanic craters. 

4) The hydrochloric acid monitoring project 
(GASPEC) moved ahead to the point the sensor is now performing 
in the laboratory. Tests will be run to determine how well 
and then further reduction of the noise will be attempted. 

This further work depends on hiring a new technician which is 
planned for early 1979. 

5) Funding for a planned expedition to the 
Antaractic to do volcano research with Ohio State University 
personnel in December 1979 was denied by the NSF. We await 
copies of the reviewers comments before deciding our next move. 
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6) We expect to visit Mt . Etna once in 1979 to 
establish base levels of output of sulfur dioxide. A paper 
relative to premonitory increases in sulfur dioxide at Etna 
volcanic eruptions has been submitted for publication to 
Nature . 

7) Work at both Kilauea and Mauna Loa Volcanoes 
in Hawaii has been stepped up. The first visit to Kilauea with 
the COSPEC was in January 1975. Very similar output of sulfur 
dioxide was recorded on our latest visit October 1978. At 
Mauna Loa in 1978 the output was at most a few tons per day. 

This visit was sponsored by this grant, the United States 
Geological Survey and the NOAA people. (See Trip Report, 
Malinconico, Stoiber September 29-October 17, 1978). Plans 
are in progress for the continuous monitoring of Kilauea for 

a month in July or August 1979 in collaboration with USGS. 
Further investigations of Mauna Loa will be instituted as part 
of our study of the major components of Hawaiian volcanic 
gases. We hope to study condensates of Hawaiian fumaroles 
during the coming summer. This will be a master's degree 
project for a Dartmouth student if present arrangements go 
forward. Since the Hawaiian volcanoes are genetically different 
from those of circumpacif ic and Kilauea is also likely to erupt 
every few years, the studies in Hawaii will be important in our 
overall program. 
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8) We were invited to give a paper relating the 
information regarding Hawaiian volcanic gases to the study 
planetary atmospheres. The information in the paper was perhaps 
reported somewhat prematurely but we took this opportunity to 

be represented in the JPL Symposium on Planetary Atmospheres 
in January 1979. 

9) The immediate goal of our gas program is to 
determine the SC> 2 output at different volcanoes and at different 
stages in volcanic activity. This requires the piecing together 
of innumerable pieces of data and will perhaps require some 
continuous monitoring for periods of a month or so at selected 
places. It is in this regard that we have chosen to monitor 
Kileauea. It appears that, at the present the state of our 
information, the output of sulfur dioxide during the various 
stages of activity, changes by the order of magnitude or more, 
where as the differences between several volcanoes at the same 
stage is decidedly less than half an order of magnitude. The 
goal, the overall pattern of volcanic emission at volcanoes of 
several kinds, will not be reached in a single year but will 
demand several years study. Our new Hawaiian work described 
above is another facet of this overall press for information 
that can be integrated into a world volcanic sulfur dioxide 
budget . 
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Small Forest Cuttings Mapped With Landsat Digital Data 
by Emily Bryant, Arthur G. Dodge, Jr., and Martha Jean Eger 


NASA Science Grant 5014 supports a cooperative project between the 
Goddard Institute for Space Studies, New York City, Cooperative 
Extension Service, University of New Hampshire, and Dartmouth College, 
Hanover, New Hampshire. Our goal is to produce useful forest information 
through computer classification of Landsat data. We have mapped forest 
types (hardwood and softwood) with some success. Forest types do not 
make up the entire picture, however. A forester also needs frequent 
updating on the location, size, and condition of harvested areas (cuttings) 
to manage his lands. We have thus turned our efforts toward mapping 
forest cuttings, specifically clearcuts in northern New Hampshire. 

Using Landsat data of June 26, 1975 (scene 5068-14433) , we developed 
a supervised classification of clearcut areas. A total of 60 cuts in 
northern New Hampshire were used as training and test areas. Cuts in 
this region range in size from 5 to 500 acres, commonly about 40 acres. 
Different stages of regrowth in the cuts had different signatures; we 
found we were able to distinguish five general age categories. 

We provided public and private foresters with sample printouts of 
their land for evaluation. They were satisfied with the shape and location 
of the cuts. Accuracy of size compared to ground truth provided by field 
foresters was variable, but to quote one, "probably your figures are 
better than ours". 

An essential part of mapping clearcuts is identifying changes from 


To C r p re 3^/c J / / > riv $y ^ j , 

' i/m.<k, /)<>'•'! /*?* : £'' cl 1 


f 


2 


v 


year to year. This requires classification categories which are consis- 
tent from one Landsat pass to another. We are currently experimenting 
with the development of consistent forest categories for three Landsat 
passes over a base area with little change during the time of the passes 
Difference in acreage estimates from the three passes in three forest 
categories was rarely more than 15% in selected 30,000 acre areas. 

We feel that the capability to map clearcuts brings Landsat a step 
closer to being an operational part of current forest inventories and 
management activities. 



A KEY FOR BETTER INTERPRETATION OF LANDSAT DIGITAL DATA 
• ■ ' by - 

Arthur G. Dodge, Jr.., Emily S. Bryant 6 Martha Jane Eger 

Tools should be developed that increase accuracy, save time and improve or 
broaden the skills of people involved in remote sensing activities. 

Keys are tool's that have been used by photo interpreters for many years. They 
can be descriptive, pictorial , or both, in format. Their usefulness is only 
limited by the imagination of the creator or the ability of the user.. 

Adapting Landsat digital data into established inventory systems or applying . 
satellite data for solving. forest management problems' has often been questioned by 
field users and potential users. 

We have observed that one road block to more use of Landsat digital data is a 
lack of understanding, by the user, of what is being displayed on a printout or 
CRT screen. Another observation is that effective interpretation of Landsat data 
increases in direct proportion to correlation of that data with conditions in the 
field. 

Field checking the interpretation or analysis of any remote sensing data is 
exacting, tedious, time consuming and expensive. It can also be fun or very un- 
comfortable. Historically, field checking time has become more limited due to 
financial restrictions and added new demands upon all users of remote sensing data. 

We have developed a photo-digital key which graphically describes-forestry field 
conditions and illustrates how these conditions are displayed on a Landsat digital 
printout. The key was developed while mapping forests in northern New Hampshire. 

We have been able to include in the key: 

Water 

Hardwood Forest 
Softwood Forest 
Mixed Wood Forest 

Clearcuts in hardwoods or mixed wood with five stages of forest regrowth 
Clearcuts in softwood forest 


•t n C I 2 -"" 



- 2 - 


We believe that this key can be used to analyze Landsat data covering northern 
New England, Slight modification could adapt it to other areas of the Northeastern 
United States and Canada. 

. Advantages of this key are: 

Introductory description of forest conditions for beginning data analysts. 
More accurate interpretation of Landsat forest mapping before field checking. 
Less field checking time needed. 

A better way to describe Landsat forestry mapping to new users and 
potential users. 


NOTE: This is a joint project between the Cooperative Extension Service. 

University of New Hampshire. Dartmouth College , and Goddard Institute 
for Space Studies, New York. MY, -under NASA Grant Sort. 
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APPLICATION OF REMOTE SENSING TO RECONNAISSANCE GEOLOGIC 
MAPPING AND MINERAL EXPLORATION 


Richard W. Birnie and Jon D. Dykstra 
Department of Earth Sciences 
Dartmouth College 
Hanover, NH 03755 


ABSTRACT 


A method of mapping geology at a reconnaissance 
scale and locating zones of possible hydrothermal 
alteration has been developed. This method is based 
on principal component analysis of LANDSAT digital 
data and is applied to the desert area of the Chagai 
Hills, Baluchistan , Pakistan. A method for airborne 
spectrometr ic detection of geobotanical anomalies 
associated with porphyry Cu-Mo mineralization at 
Heddleston, Montana has also been developed. This 
method is based on discriminants in the 0.67 urn and 
0.79 urn region of the spectrum. 


I. INTRODUCTION 

This paper considers two applications of remote sensing to geology. ■ The 
first method uses computer processed LANDSAT digital data to map geology over 
large areas of well exposed surface features. The second method employs a. 
high resolution airborne spectrometer to detect geobotanical anomalies asso- 
ciated with Cu-Mo mineralization. 

II. DIGITAL PROCESSING OF LANDSAT DATA 

A computer algorithm for the classification of LANDSAT digital data has 
been developed at Dartmouth College in cooperation with the Goddard Institute 
for Space Studies, New York. This algorithm has been applied to the Chagai 
Hills, which form an eliptical massif covering about 170 by 40 km in south- 
western Pakistan (Figure 1). The region is well suited for a remote sensing 
study because of its extreme aridity and near total lack of vegetation which 
might otherwise contaminate the spectral signatures of the surface geology. 
Furthermore, due to the rugged topography and limited access to interior 
regions of the Chagai Hills, remote sensing is a cost and time efficient way 
to gain information on the extent and distribution of surface geologic types. 

The Chagai Hills are a predominantly igneous terrain. The major forma- 
tion is the Sinjrani Volcanics which consists of over 1300 m of andesitic 
flows, agglomerates, volcanic conglomerates, and pyroclastics with small 
amounts of rhyolite and basalt (Hunting Survey, 1960). The intrusive rocks 
in the Chagai Hills, the Chagai Intrusives, are divided into two groups: a 

dark, main phase of granodiorite-quartz diorite composition and a lighter toned, 
younger phase of granitic composition. 

The Chagai Hills are believed by these authors to be an extension of the 
Tethyian Suture Zone of Iran. Several porphyry occurrences are known in Iran 
along the Suture Zone; the Saindak porphyry copper deposit lies west of the 
Chagai Hills, and new prospects in the area between the Chagai Hills and 
Saindak have been located using computer processing of LANDSAT digital data 
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(Schmidt and Bernstein, 1977). Several prospects are known to exist in the 
Chagai Hills and the possibility of additional occurrences is believed to be 
good . 


The classification algorithm differentiates seven surface geologic types 
over the 5200 m2 of the Chagai Hills. These include the bedrock units de- 
scribed above: 1) the Sinjrani Volcanics, 2) the granodiorite-quartz diorite 

phase of the Chagai Intrusives, and 3) the granitic phase of the Chagai 
Intrusives. Three unconsolidated surface units are also mapped and these 
include: 1) volcanic float derived from high standing Sinjrani Volcanics and 

often mantling low lying Chagai Intrusives, 2) dry wash deposits of clastic 
sediments, and 3) granitic talus localized aroi high stands of Chagai 
granites. The last surface geologic type is zones of hydrothermal alteration* 
which are characterized by areas of subtle iron oxide staining created by 
weathering of sulfide minerals associated with phyllic and potassic altera- 
tion. 


The classification of the digital data requires choosing control areas of 
the surface geologic types which are to be mapped. These were chosen on the 
basis of field and aerial photographic studies. The distribution of reflected 
radiances of pixels of a given control area type is represented by- a "charac- 
teristic line" in the LANDSAT 4-dimensional color space. This "characteristic 
line" is a line parallel to the principal component of the data and passing 
through the center of gravity of the data. Examples of "characteristic lines" 
for hypothetical distributions of pixels within two control area types 
(A and B) are plotted in two dimensions in Figure 2. In the actual classifi- 
cation procedure, there is a 4-dimensional "characteristic line" for each 
control area type. 

The LANDSAT color space can be sliced by lines of constant total reflec- 
ted radiance representing the sums of the reflected radiance values in the 
four LANDSAT bands. These lines of constant total reflected radiance .are 
shown as dashed lines in the 2-dimensional representation in Figure 2. These 
lines divide the color space into "slices" with slice 01 representing total 
reflected radiance values achieved only by control area type B; slice 82 in- 
cludes total reflected radiance values of control areas A and B; and slice S3 
contains only A. This slicing is equivalent to defining areas of overlap and 
non-overlap of the distribution of pixels as projected onto a line perpendicu- 
lar to the lines of constant total reflected radiance. The LANDSAT 4-dimen- 
sional color space can be sliced in a similar manner to determine which con- 
trol areas are to be considered possible classifications based on the total 
reflected radiance value of an unknown pixel. An unknown pixel need only be 
tested for possible classification with those control area types whose 
distribution of total reflected radiance values overlaps that of the unknown 
pixel . 


The. next discriminant in the classification procedure is based on an 
unknown pixel's distance from the "characteristic lines" of the individual 
control areas.. These distances are measured on the perpendicular from the 
unknown pixel to the "characteristic line" (DEL/; or DEL/j in Figure 3). The 
separation of two "characteristic lines" is also calculated (DIF/\j 3 in Figure 3 \ 
The likelihood of a given pixel's belonging to a control area type is quanti- 
fied by calculating a "confidence value" defined as DIF/DEL. The "confidence 
value" is, therefore, directly proper tional to the separation of the 
"characteristic lines" (DIF) and inversely proportional to the pixel's dis- 
tance from the "characteristic lines" (DEL). Classification information in a 
single pixel is augmented by that of its neighbors by summing "confidence 
values" of a central pixel and its eight contiguous neighbors. The central 
pixel is classified as the control area with the highest summed "confidence 
value" . 

The above described classification procedure applies to reconnaissance 
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mapping of regional geology where each control area has a sufficiently large 
pixel population to establish a well defined "characteristic line". Known 
areas of hydrothermal alteration in the Chagai Hills provide only five pixels 
which satisfactorily represent the zones of potential mineralization. A 
reliable "characteristic line" cannot be calculated from so few datum points, 
so the classification of zones of possible hydrothermal alteration is based on 
individual pixel spectra. Pixels in the Chagai Hills that have spectral 
signatures exactly matching those of the known hydrothermal alteration pixels 
are classified as altered if they and at least two of their eight contiguous 
neighbors have hydrothermally altered signatures. 

Every pixel in the Chagai Hills was classified by the above procedure, 
and a 1:48,000 scale map of surface geologic types was generated. This final 
LANDSAT geologic is too large to be included in the format of these proceed- 
ings; however, portions of this map will be shown and compared to portions of 
the 1:253,400 scale photo-geologic map of the Hunting Survey (1960). A 
portion of the Hunting photo-geologic map is shown in Figure 4a. This map 
distinguishes zones of Sinjrani Volcanics (SjV) , volcanic float (Rec) , and 
zones of Chagai Intrusive (Chi). The LANDSAT geologic map for the same area 
is shown in Figure 4b. The photo-geologic map does not differentiate the 
intrusive phases, but the satellite identifies the intrusive in this area as 
the older granodiorite-quartz diorite phase (Chi (a)). This classification 
agrees with our field studies. The agreement between the photo-geologic map 
and the LANDSAT map is seen to be good, particularly when the geologic inter- 
pretation based on our field observations is made that the area mapped as 
volcanic float (Rec) on the LANDSAT map is underlain by Chagai Intrusives. 
Another area of the Chagai Hills is shown in Figures 5a and 5b. The LANDSAT 
geologic map distinguishes two intrusive phases: 1) granite (Chi (b) ) and 

2) granodiorite (Chi (a)). The Hunting Survey's (1960) photo-geologic study 
did not separate these two phases, and showed the Chagai Intrusives as a 
single, undifferentiated unit (Chi). The dimensions of the granite shown in 
the LANDSAT map agree well with our field observations. Zones of possible 
hydrothermal alteration have been identified in regions of the Chagai Hills, 
and these zones will be the subject of upcoming field studies. 

It should be emphasized that the entire 5200 km 2 of the Chagai Hills have 
not been field checked. The lack of detailed ground truth makes estimation 
of errors of omission and comission impossible. However, we draw a largely 
favorable qualitative comparison between our LANDSAT surface geologic map and 
the photo-geologic map of the Hunting Survey (1960). It should be further 
emphasized that a good deal of geologic information which is impossible to 
obtain from the LANDSAT mapping procedure described herein can be derived from 
conventional field techniques. However, we believe that this LANDSAT geologic 
mapping procedure offers a valuable method of reconnaissance geologic mapping. 
The main advantage is that it allows reconnaissance scale mapping of large 
regions at a fraction of the expense and man hours involved with conventional 
techniques. Another advantage over conventional photo-interpretation studies 
is the ability to make geologic discriminations on the basis of color. In 
addition, the objective nature of the machine processing decreases the 
possibility of subjective errors. 

III. AIRBORNE SRECTROMETRIC DETECTION OF 
GEOBOTANICAL ANOMALIES 

Where surface exposures of the geology are masked by vegetation, an 
airborne spectrometer has been used to detect cjeohotanical anomalies. These 
anomalies manifest the toxic effects of mineralization and thus may be used 
as targets in an exploration program. The airborne spectrometer used in this 
study was developed by William Collins and is described by Collins (1976). 

The instrument measures reflected radiance (mW/cm 2 .sr) in the visible and 
near infrared region of the spectrum. The spectral band width is about 
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0.0014 pm. The instrument was mounted in a Rockwell Aero Commander 680-E and 
flown at an elevation of about 600 m above the surface. The reflected ra- 
diance was integrated over a surface area about 18 m square, and photographs 
were taken every 10 spectra to permit a correlation of the reflected radiance 
spectra with the ground surface. 

The area chosen for this study is the Heddleston copper-molybdenum 
porphyry deposit located about 50 km NW of Helena, Montana. The geology is 
reported by Miller et al. (1973). The deposit is localized in a complex 
series of tertiary quartz monzonite and quartz porphyry intrusives. These 
prophyries intrude the Precambrian Spokane Formation of the Belt Series and 
are cut by a Cretaceous (?) diorite sill. A broad 1.6 km2 area of the in- 
trusives shows sericitic, silicic, and argillic alteration. The mineraliza- 
tion includes primary chalcopyrite and molybdenite with zones of supergene 
chalcocite and covellite. Drilling has delineated zones of shallow to inter- 
mediate depth mineralization (Figure 9) , which are thinly covered by alluvium 
and leached capping (Miller et al., 1973). Dense pure "stands of lodgepole 
pine (Pinus contorta) up to 25 m high cover much of the deposit. 

The deposit was traversed by a series of flight lines that covered both 
the mineralized and nonmineralized zones. Each flight line was about 2 km 
long and included about 100 spectra, of which between 10 and 30 were of 
lodgepole pine from the mineralized zones (Figure 9). 

The first test was to determine if spectra known to be from the 
mineralized zone were different from spectra from outside the mineralized 
zone. Photos taken from the aircraft were used to isolate spectra of uniform- 
ly dense lodgepole pine. Spectra that included areas of streams, roads, grass, 
or nonuniformly dense trees were rejected. The remaining spectra were 
separated into one of four categories depending on whether they were sunlit, 
shaded, mineralized, or nonmineralized: 1) sunlit- nonmineralized, 2) sunlit- 

mineralized, 3) shaded-nonmineralized , and 4) shaded-mineralized . The spectra 
in each category were then averaged and compared. Figure 6a shows an average 
of 10 suniit - mineralized and 10 sunlit - nonmineralized spectra from flight 
line SI. The mineralized spectra have higher reflected radiance values in the 
visible portion of the spectrum. In the near infrared portion, the mineralized 
reflected radiance values are also higher than the nonmineralized values, but 
proportionally the difference is less. By dividing the mineralized by non- 
mineralized spectra channel by channel, the two average spectra are compared 
in a ratio plot (Figure 6b) . The greatest proportional difference in the 
average spectra occurs in the 0.67 pm region of the spectrum. The difference 
is least in the 0.79 urn region. Eight flight lines, both morning and after- 
noon, yielded sufficient numbers of spectra for 11 such ratio plots. Of 
these 11 ratio plots, four compared sunlit spectra and seven compared shaded 
spectra. Nine of the 11 plots showed the general properties depicted in Figure 
6, i.e., higher 0.67 um and near equal 0.79 ym reflected radiance values for 
the mineralized spectra. Two ratio plots showed equal or lower 0.67 ym 
reflected radiance values for the mineralized spectra, but equal values in the 
0.79 ym region (Birnie and Hutton, 1976). 

The experiment described above required prior knowledge of the zones of 
mineralization; therefore, it does not provide an exploration tool. A method 
of analyzing the spectra without prior knowledge of the geology has been 
developed employing data representing an average of 12 spectrometer channels in 
the 0.67 ym and 0.79 ym spectral regions. First, a flight line of data is 
analyzed spectrum by spectrum to determine which spectra represent uniformly 
dense vegetation; those spectra representing nondense vegetation such as roads, 
grass, etc., must be removed. Spectra that do not represent uniformly dense 
vegetation have anomalously high reflected radiance values in the 0.67 ym 
region. Figure 7a shows a plot of the 0.67 ym reflected radiance values across 
flight line #11. Those values that are associated with streams, roads, grass, 
and nondense vegetation are indicated by a cross. All but one of these 
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"bad spectra" have reflected radiance values greater than 0.00125 mW/(cm 2 .sr. 
channel). It is concluded, therefore, that it is possible to discriminate 
between dense vegetation spectra and other . spectra by rejecting all. spectra 
that have 0.67 pm reflected radiance values greater than a particular cutoff, 
in this case 0.00125 mW/ (cm^ . sr . channel ) . 

The dense tree spectra are separated into sunlit and shaded categories 
based on a discriminant at 0.79 yin where the effects of shading are emphasized 
by a relatively strong reduction in reflected radiance. It should be recalled 
that the effects of mineralization on the spectra are minimal in the 0.79 u m 
band. In Figure 8, the 0.79 |jm reflected radiance values are plotted for 
flight line #11 with the sunlit and shaded zones indicated. If a cutoff is 
applied at the 0.00375 mW/ (cm 2 . sr. channel) and all spectra with reflected 
radiance values above and below this cutoff are considered sunlit and shaded, 
respectively, a good correlation with the ground truth is achieved. 

Now it is necessary to discriminate between the mineralized and non- 
mineralized spectra within each of the shaded and sunlit categories. The 
ratio plot in the averaging experiment (Figure 6) indicates that the 
mineralization effects are most pronounced at 0.67 pm. Figures 7b and 7c 
show the 0.67 |im plots of reflected radiance values for the sunlit and shaded 
spectra, respectively. A cutoff can be applied to these data with all values 
greater than 0.0085 mW/ (cm2 . sr . channel ) i n (-he sunlit zone and 0. 0070 mW/ 

(cm 2 . sr. channel ) in the shaded zone classified as mineralized. In this 
example, most of the mineralized spectra are from a sunlit region; however, 
the same line was flown in the afternoon when the mineralized spectra were 
shaded, and a similar cutoff procedure isolated the mineralized zone within 
the shaded spectra. 

If the zones of mineralized and nonmineral ized spectra are recombined as 
in Figure 7d, a comparison can be made between the zone of mineralization as 
identified by the anomalous 0.67 ym reflected radiance values and the zone of 
mineralization identified by Miller et al. (1973). The correlation is good 
with two errors of omission within the mineralized zone and several errors of 
comission in the nonmineralized zone. 


When the eight flight lines of data that traverse the Heddleston deposit 
are analyzed in a similar manner, the results plotted (Figure 9), and the zone 
of anomalous 0.67 ym spectra enclosed, it is seen that the anomalous spectra 
identify a zone that includes the shallow to intermediate depth mineralization 
but also extends beyond this zone. It is reasonable to expect that vegetation 
would manifest anomalous reflectance values beyond the zone of mineralization, 
because any soil geochemical anomalies would be expected to be dispersed beyond 
the border of the bedrock mineralization. The errors placing mineralized 
spectra outside the mineralized zone are particularly noticeable in the 
eastern part of the map. These erroneous spectra appear in areas where the 
vegetation is anomalously sparse. 


While this method needs testing both 
other types of vegetation, the results of 
spectrometer measuring reflected radiance 
potentially powerful exploration tool for 


in other mineralized areas and with 
this study indicate that an airborne 
values of vegetation represents a 
the economic geologist. 


799 


ACKNOWLEDGEMENTS 


Robert Schmidt of the U.S.G.S. provided invaluable assistance in the 
LANDSAT and Pakistan portion of the project, and Dr. Victor E. McGee of 
Dartmouth College provided extensive assistance with the principal component 
analysis. The field work in Pakistan was supported by National Science 
Foundation grant “INT 77-07767 and by the Resource Development Corporation 
(RDC) of Pakistan. The help of RDC's general manager, Dr. S.A. Bilgrami, is 
greatly appreciated. Dr. William Collins of Columbia University and the 
Goddard Institute for Space Studies (GISS) , New York, developed the airborne 
spectrometer. Dr. Stephen Unger and his staff at GISS, along with Dr. Gerry 
Carlson of the Department of Earth Sciences, Dartmouth College, provided 
extensive help with the data reduction. As the primary research assistant on 
the Heddleston project, Michael Hutton of the Department of Earth Sciences, 
Dartmouth College, provided particularly important contributions to the 
computer analysis of the spectral data. The computer processing of the data 
and the field work at Ueddleston, Montana were supported by NASA grant fiNSG 
5014. The Anaconda Company allowed us access to their geologic information 
on Heddleston. 


REFERENCES 


Birnie, R.W. and M.S. Hutton, 1976, Reflective spectra of vegetation growing 
over the Heddleston copper-molybdenum deposit, Montana (abs.): Abstracts 

•with Programs, G.S.A. An. Meet., Nov. 1976, Denver, p. 779-780. 

Collins, W. , 1976, Spectroradiometr ic detection and mapping of areas enriched 
in ferric iron minerals using airborne and orbiting instruments: Ph.D. 

Dissertation, Columbia University. 

Hunting Survey Corporation, Ltd., 1960, Reconnaissance Geology of Part of 
West Pakistan: a Colombo Plan Cooperative Project: Toronto, 550 p. 

Miller, R.N., E.P. Shea, C.C. Goddard, Jr., C.W. Potter, and G.B. Brox, 1973, 
Geology of the Heddleston copper-molybdenum deposit, Lewis and Clark County, 
Montana: presented at Pacific Northwest Metals and Minerals Conference, 

Coeur d'Alene, Idaho, April 12-14, 1973. 

Schmidt, R.G. and R. Bernstein, 1977, Evaluation of improved digital pro- 
cessing techniques of LANDSAT data for sulfide mineral prospecting: in 

Proc. First An. William T. Pecora Mem. Symp . , Oct. 1975, Sioux Falls, 

U.S. Geol. Surv. Prof. Paper 1015, p. 201-212. 


800 



LANDSATV 
MOSAIC 


100 KM 


Figure 1. LANDSAT MOSAIC 

Mosaic of six LANDSAT frames covering the study area. The 
study area is that portion of Chagai Hills which is in 
Pakistan (cross-hatched section of the eliptical feature 
in inset location map) . 




Figure 2. CHARACTERISTIC LINES 

Distribution of pixels of two 
hypothetical control areas 
showing "characteristic lines" 
(A+B) and lines of constant 
total reflected radiance (dashed) . 


Figure 3. DIF AND DEL 

An unknown pixel's (dot) 
perpendicular distance from 
each control area (DELa and 
DELg) and the separation of 
the two control areas (DIFab) 
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Figure 9. MAP OF MINERALIZATION AND ANOMALOUS SPECTRA 


The zones of shallow to intermediate depth mineralization 
at the lleddleston deposit are indicated by shading (after 
Miller et al., 1973). Traces of individual flight lines 
arc shown with snectra classified as mineralized and 
nunninera l L/.ed indi cated hv open triangles anil dots 
respectively. The solid bold line interprets the 
mineralized zone from the spectral data. 
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ABSTRACTS WITH PROGRAMS. 1978 



IV, and hy local contact Ion exchange In VI , suggesting presence of per- 
vasive fluid In Jb-IV, but absence in VI, conaintent with features of 
cataclasltes in these zones. 


waters. 7) Oxidation of iron in wall rocks may be an effect of, rather 
than a cause for, uranium precipitation. 8) Pitchblende can precipltatf 
under conditions where pyrite or hematite are stable, or where iron Is 
mobile, resulting in the bleaching of wall rocks. 


LONG-O I STANCE CHRONOCORRELATION IN THE MIDDLE CAMBRIAN 

Robison, R. A., Department of Geology, University of Kansas , 
Lawrence, Kansas 66045; McGee, J. W. , Cities Service Oil Company, 

Box 300, Tulsa, Oklahoma 74102; Resell, A. J., Department of 
Geology, University of Kansas, Lawrence, Kansas 66045; Strickland, 

0. K. , Department of Geology, University of Kansas, Lawrence, Kansas 
66045; and White, W. W. , III, U.S. Bureau of Land Management , 

Buffalo Resource Area, Buffalo, Wyoming 82834. 

Most Cambrian species are endemic, which creates severe problems in 
global chronocor re I at Ion . The majority of agnostoid trilobltes were 
probably pelagic and some were cosmopolitan. Two species, 

Ptychagnos tus gi bbus and P. atavus , show widespread stasis throughout 
their geographical range, but both are separated by morphological gaps 
from their putative ancestors. Their phylogeny is consistent with, a 
punctuated-equi 1 I br ium model of speclation. 

The bases of both the P. gibbus and P. atavus zones have previously 
been suggested as potential horizons for series boundaries. In many 
parts of the world the Incoming of P. g i bbus is associated with a 
marked change In lithology, and In detail, t^e initial appearance of 
the species may have been delayed at any given locality. In contrast, 
P. atavus commonly occurs in monofaclal open-ocean successions and 
probably Its Initial appearance closely followed speciation. In 
western Utah, a likely area for establ 1 shment of a series boundary- 
stratotype, P. g i bbus and other trllobites occur with inarticulate 
brar.hlopods Tn the" lower part of the Wheeler Shale, which accumulated 
under relatively deep-water conditions'. Without significant change in 
either lithology or the inarticulate brachiopod assemblage, P. g i bbus 
is abruptly replaced by P. atavus. Associated changes In poTymeroid 
trllobites, although less extens i ve , resemble changes at the bases of 
Late Cambrian biomeres ..which have been Interpreted as adjustments to 
cooler water. Considering all factors, chronocor re I at i on based on the 
first appearance of P. atavus Is probably as precise as any in the 
Cambr i an . 

TRACE ELEMENT GEOCHEMISTRY OF MAFIC AND FELSIC NINETTES , BUELL PARK 
HATREME, NAVAJO VOLCANIC FIELD 

RODEN, Michael F., Department of Earth and Planetary Sciences, 
Massachusetts Institute of Technology. Cambridge, Massachusetts 
02139, and Department of Ceologlcal Sclencea, University of Texas, 
Austin, Texas 78712 

Potass lum-r Ich nineties (3.5-7Z K 2 0. K 2 0>Na20) of the Navajo volcanic 
field containing abundant spinel and rarer garnet lherzolite Inclu- 
sions were erupted 25-30 million years ago in the central Colorado 
Plateau. Samples of mafic (49Z Si0 2 , 10Z MgO) to felslc (60Z S10 2 . 

6Z MgO) mlnette from the Buell Park composite klmbertitlc breccla- 
ralnette diatreme were analyzed for rare earth elements (REE), Rb, Ba, 
Sc, Cr , and Ta by instrumental neutron activation. Chondrite- 
normal ized REE patterns for all the minettes are similar and highly 
fractionated (La - 330-370 x chondrlte, Yb - 5-7 x chondrite). Also, 

Ba Increnaps slightly through the suite (2100-2350 ppm), while Rb 
increases by two-fold (100-200 ppm). In contrast, compatible clemcnta 
decrease significantly from mafic to fclaic mlnette: Cr, 410 to 100 

ppm, and Sc, 16 to 9 ppm. A major element model based on whole rock 
and phenocryat analyses relates mafic to felsic mlnette by fraction- 
ation of dlopslde, phlogopite, magnetite, and apatite. Apntlte frac- 
tionation was important in controlling the REE and explains decreasing 
P 2 05 and P 2 05 /Ce through the aulte. The fractionation of Cr and the 
"buffering" of the REE are not easily explained by a model in which 
mafic and felsic mlnette are derived by different amounts of partial 
melting of the same source. P 2 05 /Ce ratios (^ 40) of mafic mlnettca 
are lover than that suggested for the source of alkali basalts and 
nephcllnltes and are consistent with a source, for the minettes 
enriched in light REE. A garnet lherzolite which initially contained 
apAtlte and phlogopite la a possible source. 


HYDROTHERMAL TRANSPORT AND DEPOSITION OF URANIUM, AND THE ORIGIN OF 
VEIN URANIUM DEPOSITS. 

ROMBERGER, Samuel B.. Department of Geology, Colorado School of 

Mines, Golden, Colorado 80401. 

The stabilities of various aqueous uranium complexes are evaluated at 
up to 300 l using available thermodynamic data. The distribution of 
aqueous uranium species are plotted as a function of pH and log oxygen 
fugacity to determine the complexes responsible for the transport, and 
the mechanisms for deposition, of uranium in natural hydrothermal 
systems. When these results are considered along with the hydrothermal 
alteration and mineral assemblages observed in vein uranium deposits, 
the following conclusions are reached: 1) Carbonate complexes become 
less important at elevated temperatures and are not significant in the 
transport of uranium in hydrothermal systems. 2) In deposits containing 
fluorine, uranyl fluoride complexes are most important in uranium 
transport. 3) Chloride complexing is not significant even where 
chloride activites exceed 1. 4) Under certain conditions phosphate 
complexing may be important in transporting uranium. 5) The solubility 
of pitchblende decreases as temperature increases, and therefore 
cooling is not an adequate mechanism for deposition. 6) Uranium is 
deposited by an increase in pH or a decrease in activity of complexing 
anions. These changes may occur as a result of reaction with wall 
rocks, boiling and loss of volatile components, or mixing with ground 


LARGE-SCALE "GRENVILLE" NAPPE STRUCTURES INVOLVING GRANITIC GNEISSES If 
SOUTHEASTERN ONTARIO AND THE NORTHWEST ADIRONDACK LOWLANDS 

ROME'Y , Dill and JACOBY, Russell S., Department of Geology and 
Geography, St. Lawrence University, Canton, H.Y. 13617 
A competent, rumpled, northwest-dipping sheet of granitic gneiss 
extends from Parishvllle to Alexandria Bay, N.Y. Apparently Isolated 
leucogneisses (the Canton, -Gouverneur, Hyde School, Fish Creek, 
California, Clark Pond, and other "phacol 1 thic" bodies) result from 
"egg-carton" structure of the sheet. Contorted, plastically deformed, 
incompetent metasedimentary rocks fill troughs and depressions between 
domal arial culminations on the granitic sheet. 

Structural relationships show the Parishville body, to the east, 
as the hinge of a major, recumbent, isoclinal nappe, with one limb 
extending into the Adirondack Highlands. To the west, the granitic 
gneisses at North Harrmond and Alexandria Bay extend the sheet across 
the St. Lawrence River to include the Rockport granite in Ontario, 
possibly forming another major fold hinge. The amplitude of these 
folds approaches 60 miles. They. extend a competent granitic strati- 
graphic unit across the southwestern part of the Grenville Province. 

The northeast-trending "grain" of Grenville structures represents 
axial-surface foliations to these nappes. 

Granitic rocks of the sheet may represent metasediments or 
metavolcanics of the Grenville Supergroup which have locally undergone 
melting, and remobilization. Other acidic gneisses beyond the ends of 
the fold hinges, such as charnockitic gneisses of the Diana and 
Santa Clara complexes, Now York, may be correlative with the Westport, 
Crow Lake, Gananoque, Lyndhurst. and Perth Road motaplutonlc bodies 
in Ontario and therefore represent rocks derived from anatexis of 
pre-Grenville basement rocks. 

EARLY CEN07.0IC r 1 .AT F PrOF.fi All I 7. AT I ON 

KONA, Peter A., cuul RTCMAKDGON! F.v.nn F* . , Nation.il Ocean- 
ic and Atmor.phnrio Administration, Atlant ic Oceanograph- 
ic and Met corn loF.ic.nl Laboratories, Miami, Florida 
3 3LM9 

'Also University of Miami, Roscnr.t i ** 1 f.chool of Marine 
and Atmospheric Science*, Miami, Florida, 33 lug 

Synthesis of . re*g Lorta 1 plat*' movement:; reveals a major reor- 
ganization of global plate motion from M**r.ozoic to Cenozoic 
patterns that occurred primarily dnrinr the For on** epoch 
(S3.S to 37.5 my BP). Th** reorganization involved reorien- 
tation of relative plate motions with large N-S components 
into large E-W components, cent i nua fion -t pre-existing E-W 
sea floor spreading, initiation of n*'w E-W *jpr 1 inp. , 
dece 1 ^rat ion el •;pr*'ad i nr, raters, ami obdorti-n ~ f ephio- 
litcs. The reorgan i r.at i on is at tri b*r. <*r] te increase in 
resistance impose*! on the global plat*’ system by increase 
in length of L-W trending cnlljsional plat*' boundaries fron 
7,500 to 78,509 km *luring the interval 56 to U0 my BP, 
forcing the system to reorient along lin**s <**f l**ss resist- 
ance. The pattern of global plat*? motion 1 '.-suiting from 
the reorganization is corn; * ra i nod by the resistance imposed 
by 19,000 km of colli sional plate boundaries. 

BIOTURBATE TEXTURES IN H1JDDY SEDIMENTS 

R0NAN, Thomas F.. Jr., Department of F.arth 4 Space Sciences, 
University of California, Los Angeles, California 9H026; 

L1PPS, J. H., Department of Geology, University of California, 

Davis, California 91567 

The textural properties of muddy sediments are Influenced to n large 
degree hy the mot 111 tv modes and feeding tvpes of thetr resident biota. 

A muddy sediment which contains a dense assemblage of suspension feed- 
ing Infauna which build tubes will have n Moturbate texture that 
arises primarily from the activitv of Its between- the- t ubes cemmunity. 
Myriad small sinuous burrows, primarily produced bv annelids and 
crustaceans, will pervade the sediment. Because the tubes constitute 
a barrier to subsurface locomotion, the larRc resident skeletonized 
Infauna (bivalves and shrimp) will consist of relatively sessile 
species that build permanent burrows. The overall effect Is to pro- 
duce a sediment that Is cohesive and contains well-defined biogenic 
structures. Because cemented tubes tend to degrade quick lv, pflleon- 
tnloglc recognition of 11 t iibc-hiil 1 dr r association mnv depend entirely 
upon the distinctiveness of the suite of sediment traces produced by 
the coexisting Infauna. 

By comparison, a less densely populated muddy sediment will have * 
bioturbnte texture that results from deposit feeding and locomotion 
traces. The active probing of sediment with tentacles, proboscis, 
and palps will produce a mottled to homogenous sediment texture. Th* 
overall water content of this sediment will be high and deformation 
of biogenic structures would be expected during preservation. 

1978 VOLCANIC ERUPTION CLOUD SAMPLING PROJECT 

ROSE. W. I., Jr., CADI.F. , R. D. , I.A7.RUS, A. I.., HF.IDT, L. F.., and 
D. A. GILLETTE, National Center for Atmospheric Research, Boulder, O 
80307; HUEBERT, B. J., Colorado College., Colorado Springs, CO 8090J; 
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STOIBER, R. E. , SELP, $• and BRATTON, C. , Dartmouth College, Hanover, 
NH 03755; CHUAN, R. L., Brunswick Corporation, Costa Mesa. CA 92626; 
WOODS, D. C., NASA Langley Research Center, Hampton, VA 23665; 
FRIEDMAN, l., ZIELINSKI, R. A. and SMITH, D. , U.S. Geological Survey, 
Denver, CO 80225; WILSON, L. , University of Rhode Island, Kingston, 

R I 02881. 

During February and March 1978, a team of scientists and a specially 
equipped Queen Air Aircraft flew eleven sampling flights over the 
active volcanoes of Pacaya, Fuego and Santlagulto in Cuatemnla. Through- 
out this period each of the volcanoes showed moderate activity and had 
plumes which could be readily sampled at various altitudes and distances 
from the active craters. 

Equipment on the aircraft included: 1) a treated filter sampler for 

determination of relative amounts of SOj , 1ICI, HRr, HF, SO*", particu- 
lates and trace elements; 2) a 10-stage Piezoelectric Particle Cascade 
for particulate aize distribution and alze sampling of particulates in 
the range of 25 to 0.05 iim; 3) evacuated stainless steel canisters for 
quantitative gas chromatographic measurements of HjS, COS, CO?, CO and 
SO*; 4) a SOj correlation spectrometer for the determination of SO? 
fluxes In the plumes by traversing; 5) an automatic 35 ram multi-speed 
sequence camera for dynamic photography of erupting plumes; 6) a flow- 
thru stainless steel gas sampler for H, 0 and C isotopic measurements 
of gases; 7) an IncrtlaL Navigation System for ground and wind speed 
determinations; 8) a gas bubbler sampler for trace element sampling. 

At abstract time, few Jesuits are available. 

The project demonstrates clearly the feasibility of planned missions 
to sample more or less continuous low level emission from the same 
explosive volcanoes that project material into the stratosphere in short- 
lived sporadic events. 


SNOWBALL MUSCOVITE IN THE CENTRAL SWISS ALPS: AN INTERNAL 
RECORDER OF TECTONOMETA MORPH ISM 

Rosenfcld, J. L. , Department of Earth & Space Sciences. 

University of California, Los Angeles, California 90024 
Snowball muscovite (SM) occurs with the well-known snowball 
garnet (SG) in north-dipping, s taurollte-zone , Llasslc cal- 
careous schists in the Isoclinal Plora Syncline separating 
the Gotthard and Lucomagno massifs. Its lower symmetry cau- 
ses SM to reveal petrogenetlc features beyond those shown 
by SG. Some muscovite recrystallized around new nuclei of 
scattered orientation while other coexisting muscovite re- 
acted to form snowball blotlte. Rotation of SM and their 
simultaneous, orientation-dependent, anisotropic growth- 
rate caused changes in cross-sectional shapes normal to 
their rotational axe3 from relatively equant [ ( 001 ) ~4-fol 1- 
atlon (S)]to elongate || ( 001 ) ( ( 001 ) * || S ) . Low angular velo- 
city in the latter 3tate kept S parallel to the synclinal 
axial surface (cf. D. Willis). Garnet pseudomorphs of mus- 
covite, including earlier SM, demonstrate muscovite trans- 
location by cation exchange, probably along Intergranular 
boundaries Influenced by differences In surface energy and 
changing phase equilibrium. SM, SG , and other snowball min- 
erals show rotational sense that Is Invariant with position 
In the syncline. Their south-s lde-up sense Indicates growth 
during Chadwick's tectonic Phase V after folding of the Pl- 
ora Syncline. That shear, with elongation parallel to the 
dip (e.g. boudlnage of zolsite during growth), probably de- 
rived from Ollgocene southward tilting of the bordering, 
more competent massifs, caused by southward slumping of ma- 
terial at a higher crustal level. That movement also ac- 
counts for refolding of nappe-related isoclinal folds in 
the Trlasslc dolomite at nearby Passo Campolungo. 

E-COLOGY OF LARGER, SHALLOW-WATER TROPICAL FORAMINIFERIDA 

ROSS, Charles A-, Department of Geology, Western Washington 
University, Bellingham, Washington 98225 
fodern coral reefs are inhabited by a great abundance and diversity of 
species of Foramini ferida . On exposed reef flats and exposed shallow 
intra-reef basins species of Calcarinldae and Amphi steginidae dominate 
the Thalassia (turtle grass) blades and may form the substrate of the 
basins, the upper surfaces of which are bound together by foramini feral 
pseudopods. Although most of these species are relatively small, the 
families Soritidae and Alveolinidae have representatives that reach 
shell volumes greater than 100 rim. The soritid Margi nop ora ve rtebr a l is 
occurs abundantly In a number of different reef envi ronments , includ- 
ing clinging to T halassia blades In 1 to 2 m deep protected lagoons. 

H. ve rteb ral is and other larger foramini ferids start to appear consis- 
tently on substrates composed of smaller foramini ferids in deeper 
Intra-reef basins 8 to 25 m deep on exposed reefs, and inter-reef 
channels about 25 m deep. Growth rate decreases but the final size 
necessary to reach reproductive maturity increases below a depth of 
10 m, so that the largest specimens are from the deeper locations. 

The shallow depth limits in the distribution of large reef- 
dwelling foraminlferids are related to exposure to waves and currents 
that are of sufficient strength to displace the foramini ferids and 
redeposit them In deeper water. The lower depth limits, at least for 
the soritids, are related to the depth at which their photosynthetic 
zooxanthel lae partners are no longer effective symbionts. 


Mesothelioma, a cancer of the pleura and peritoneum, occurs, with few 
exceptions, only in asbestos "trades” workers (ship insulation, tex- 
tiles, construction, etc.) and in South African crocidolite miners. 

Lung cancer is seldom found in non-smoking asbestos "trades” workers. 
The percent mortality for lung cancer and mesothelioma is plotted for 
11 groups of "trades” workers. Regression analysis of these data shows 
a linear relationship_between lung cancer and mesothelioma incidence; 
79% of the variance ( r 2 ) of lung cancer is accounted for by mesotheli- 
oma. A large part of the remaining 21% variance can be accounted foe 
by intergroup differences in smoking statistics. Cancer' mortal 1 ties 
are also plotted for adult male populations of 5 nations and for 6 
"asbestos" mininq/factory populations not exposed to crocidolite. The 
very low incidence of mesothelioma in these groups shows that they do 
not have the cancer risk experienced by those in the "trades") their 

lung cancer incidence is generally accounted for by smokin g h abits 

alone. Cancer risk to those In the ••••»( 

"trades" is attributed to exposure to both »* | | 

smoke and crocidolite. There is no evl- i ' 

dence that the mineral dust levels pre- I f „■**’ 

sently maintained in responsibly operated j 1 * 1 

mines of North America cause increased p, r . . 

Incidence of lunq cancer or mesothelioma, ' . - J t _,_. 

including those mines that process rock *‘t'' - «.».•»,»— — 

containing the commercial forms of asbes- <f * ZZL, 

tos: anthophyllite , tremolite, actinolite, j 

cummingtonite , grunerite, and chrysotile. • ' * , * . ' • * 


THE DFE PROGRAM IN ENVIRONMENTAL IMPACT OF THE DISPOSAL OF HIGH LEVEL 
RADIOACTIVE WASTE IN GEOLOGICAL FORMATIONS 

Rothschild, H. C. and Barraud, C. , Nuclear Programs Division, 
Department of Fisheries & Environment, Ottawa, Ontario, K1A 1C8. 
Fisheries & Environment Canada's interest in the problem of disposal 
of high level nuclear wastes In geologic formations involves the 
participation of two services: Envi ronmental Management Services 

(EMS) and Environmental Protection Services (EPS). EMS Is supporting 
applied research in hydrogeology and groundwater geochemistry for 
determining the transport/retardati on mechanisms of radionuclides 
in the sub-surface environment. The EPS involvement consists of 
preparing a document providing recommendations for long-term 
protection of our environment with respect to disposal of high 
level radioactive waste. This document will cover two main topics: 

(a) evaluation of long-term physical and chemical stability of the 
different rock types considered to be potentially suitable as 
high-level waste disposal media. 

(b) evaluation of the efficiency of the geologic barrier for long- 
term containment of long-lived radi onuclides . 

Conclusions drawn from this document will provide a base for 
recommending improvement of solidification techniques to produce 
synthetic matrices with long-term physical and chemical stability and 
with lover leaching rates than the present day glass or ceramic 
matrices. 


DEVELOPMENTAL CYCLE AND SYMBIOSIS OF THE RECENT NUMMULI- 
TID, HETEROSTEGINA PEFRSSSA, AND THEIR IMPLICATIONS FOR 
THE UCCOmkEWCE Ur' THE Ml CRUSPHER1C GENERATION 

R0TTGER, Rudolf, I ns ti tut fur Allgemeine Mikrobiologie 
der Universitiit Kiel, Olshausens tr . 40-60, 

23 OO Kiel, Fed. Rep. of Germany 
The developmental cycle of Hetero.s t egina denressa has been 
investigated in laboratory cuTEuren and" Tn specimens from 
natural habitats in Hawaii. As In other foram ini feran spe- 
cies, alternation of generations also exists. However, ne- 
galospheric specimens (ganont, A- form) are much more abun- 
dant than microspheric specimens (agamont, P-form). This 
seems to be true for all larger foramini fera (Recent and 
fossil). The reason for this is that multiple fission is 
the main mode of reproduction of the megalospheric gene- 
ration. Gamete formation, which leads to the microspheric 
generation, is greatly reduced. The symbiotic diatoms 
which are essential for the metaholism of this foraninifer 
are easily transferred over many negal ospheric generations 
during multiple fission. Garnets, which are rarely formed 
by negalospheric individuals, do not transfer symbionts 
to the microspheric generation. Consequently, the micro- 
spheric generation has to engulf free living diatoms 
during its early individual development and transform 
them into symbionts. The transforms t ion process makes the 
establishment of symbiosis more difficult in the micro- 
spheric generation. Such difficulties along with the 
scarcity of gamete formation may account for the small 
number of microspheric specimens found in natural habi- 
tats and for their special depth distribution. Microsphe- 
ric specimens occur to a depth of 100 m with the highest 
frequency in the deeper regions of this range. 


ASBESTOS HEALTH RISKS TO THE MINING COMMUNITIES OF NORTH AMERICA 
FOSS, Malcolm, u.S. Geological Survey, 959 National Center, 
Reston. Virginia 22092 


THE PAST AS A KEY TO THE PRESENT: PLANKTONIC FORAMINI FERA AND PTEROPODS 
IN SOUTH CHINA AND JAVA SEA SEDIMENTS 

RCGTMAN, Marcia L. , Department of Geology, Humboldt State Univ., 


Volcanology, Geochemistry, and Petrology 


panted by lava flows am) rxplrilvr vulr ; ,nt,in 
iul Ivlty. I'«irlnn our sampling, Taisya w.i^ In A 
M,itr of Intense vapor mission which begun In 
!*■*<■ t'wbrr 197? and root toned at Irani through May 
1 9 78. I. It tie or no ash van being erupted, am! 
no lova Hnw nrtlvtty vns running cither. Vhnt 
vapor rslM.iftin abated durlnt; our sampling Mights, 
magma could b<‘ directly observed through gashes 
within the Mi’Kenoy crater. hnjje lias depatled 
I row It* 70-year pattern of short -lived, frequent, 
violent pyroclastic erupt tons. Beginning tn Sep- 
tember ID? anil rmil (milnj until at least June 
19/8, Furgo has been In constant Inv-lrvel pyro- 
clastic activity. In February 197B, lava flows, 
the first since at least 1917 at Fuego, were flow- 
ing down the SW flnnk of thr rone. larger lava 
flows were seen In the E flanks In June. Mild 
Vulcan Ian activity producing substantial amounts 
id ash was Interspersed with vaporous, a«U-p»»or 
pluP*-s In February. Roth types of eruption clouds 
were sampled. ^;i. , VJ. ,l li u .IJ. 1 ! n Fe|e.m dome , which 
has been cent I minus I y erupting since I1?7. Since 
April I1/} the volcano has been producing several 
pyroclastic erupt Ions each day, some reaching 
heights •! V«. Ibirlng IVhronv this activity rnn- 
llnued In produce large, ash-rlrh eruption clouds, 
t'cntinual, relatively low level activity at rver- 
c hanging subwets of lVfttr.il America volcanoes 
allow planned missions to sample explosive vnl- 
rnno rlotttln at small ('4000 n) volcanoes In good 
(dry season) weather conditions. 

*NGAR Is sponsored by the National Science 
Founds t Ion. 
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sn.n** .I'll) 1tM.nr.FN tlONSTl tvf.nts OF VOLCANIC 
KMT I ION CUU-’tlS 

A. J. t _ J.. t *r«fr 
R. J, ihiebert 

R. 0. <adlo (nil nt: National Center for Atmos- 
pheric Research. Boulder, 00 80107*) 

V. I. Rose, Jr. (Michigan Technological Univer- 
sity. Houghton, Ml 49911 

Vi»J can I r plumes In Guntem.il;> and AJ.iwVn were 
sampled for SOj, MCI and I1F vapors «* well ns 
SO* , Cl* and l'“ In particles, Tlie mdnr ratios 
lor St's /HCI varied from 1,2 tn 7.5. and that nf 
Iin to l|F from S to 71. Selective removal of 
HCI and II F with respect to S0j as a funet Inn nf 
distance from the crater was not observed. The 
average percentages of sulfur, chlorine and 
fluorine in the particle phase are respect l ve 1 v 
2.4. IB and /.?. 

flic amount of sulfate present In the northern 
hemisphere stratosphere In Spring )97*» result- 
ing from the Fuer.o eruption In Fall, 19?/., has 
been i*s tint led from isoplelh naps of sulfate 
concent rat Ion obtained by means of WR*>7f air- 
craft and ha 1 loons. 
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AIPRORNK CORRELATION SPECTROMETER 
MEASUREMENTS OF SOj IN ERUPTION CLOUDS 
FROM GUATEMALAN VOLCANOES 

Ric h A rd E. Stoihcr 

G. DTAtton (both nt: Dopt . of Earth 

Sciences, Dartmouth Coil., Hanover, 

Nil 03755) • *-* * ‘ 

McAsuroments. of SOj emitter) during 
mild eruptions of tho Guatemalan volca- 
noes Futico and Sant i.iuu i to emitted 
in a stronn v.ijxar cloud at Pa cay a wore 
nutt)*' with an airborne remote sens i no 
correlation soectromctcr . Emission of 
300 to 1500 metric tons per day of SO 2 
is comparable with data from u round- 
based measurements made at other times 
and .at other volcanoes elsewhere. The 
airborne mode eliminates the necessity 
of roads being positioned appropr iat el y 
relative to the plume; anti allows many 
t raver Res’ to he made beneath the plume 
in a short period of time?. These 
traverses may be at any desired distance 
below the plume or downwind from the 
crater- Thus one can investigate the 
change in 50 j in the plume downwind as 
well as the instrument res)>onse at 
various distances beneath the horizon- 
tal plume due to vortical distribution 
of the gas and to atmospheric attenu- 
ation. Measurements indicate attenu- 
ation loss of at least 20? in the signal 
from a plume where it is 1.8 km above 
the aircraft and a marked drop of up to 


40> in 5i', cbntonl downwind at dint.inr«*s 
of up \tt V' km from tho volcanic era tor. 
Us* 1 *’f t it** spectrometer in an aircraft 
regurirs appropriate navigational aids 
and an exit port for the telescope with 
no overhead obstructions. 
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GAS AND H ISOTOPIC ANA USES OF VOLCANIC 
ERUPTION CLOUDS IN GUATEMALA 

W. 1. R ote. Jr. 

R. D. C.idlr 
t. F.. Heidi 
R . J .• Ihiebei t 

A. L.irrut (nil nt : National Center for Atm*<s- 
phorir Rograrrh, Boulder, CO R01P?*: pte*. 
address (Rose): Michigan 1 ec hno log 1 r j 1 Univer- 

sity, Untight on. Ml 49911 
I. Friediun (U. S. Ceologlcal Survey, Denver CO 
80275) 

Cna samples were collerted by slrcraft entering 
volcanic eruption clouds of three Gu. 1 teaBl. 1 n 
volcanoes. Gas Chrooi t ogr *ph I c analyse* show 
higher H. and S gjg content* in a*h eruption 
clouds and lover H^ and S gases In vaporous gas 
plumes. H 1 60 top Ic data demonstrate lighter 
Isotopic distribution In water vapor of ash 
eruption clouds than In vaporous gas pluses. 

H***t of the H^O in the vaporous pluses Is prob- 
ably oeteoric. The data tre the firm direct 
gas analyses of explosive eruption clouds, and 
demons! rate that. In spite of atmnsphr r 1 c admix- 
ture, useful compos 1 1 1 onal information on erup- 
tive gnsrs tnn be obtained using aircraft. 
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PARTICLE SIZE DISTRIBUTION AND 
VORPIIOLOGY OF ASH IN THE FUERO AND 
SANT I AGU 1 TO PLUMES 

Stephen Self (Pcpt- of Earth Sciences, 

Dartmouth College, Hanover, NH 0)755) 

Airborne collection of pyroclastic 
p.ir t icul.it os in the plume of Fuego and 
Snntiogu ito Volcanoes, Guatemala was 
conducted in February 1978. Particles 
were collected on both Fluropore filters 
and on a gunrtz crystal microba lance 
(yCM) cascade impact or between 12 and 
15,000 ft. above r.en level am) at 
various distances downwind. SKM exami- 
nation of particles shows size ramie, 
particle type and morphology. Solid 
particles in tho plume ore juvenile 
scoria, crystal fragments anti comminuted 
non-juvenile lithic material. Sizes 
range from 1mm to <0.05ui» and the prog- 
ressive downwind decrease of the coarse 
fragment content can be documented. 

This is compared to plume dynamics obtai- 
ned from movie film of the explosions 
which produced thr ash contained in the 
plume. Initial plume velocities of 
‘10 m/s, which decelerate dramatically 
to a few meters per second before being 
blown downwind, supported particles 
<25 ura up to distances of *15 km from 
Fucqo Volcano with bulk concentrations 
of several thousand vg/m*. Data wi 1 1 - 
also be presented on the relative concen- 
tration of the various particle types 
and the total grain size distribution 
in an eruption plume produced by mild 
vulcanian explosions. 
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SIZE 0IS1R1RUTI0N, MMRriltH.nCY AND CMMKNIAI. 
COMPOSITION OK FINE PARTICLES IN VOICANIC PLUMES 

|). •!. V.»,hIm {NASA Lmr.lrv Research Center, 
!l;m|<t«n, VA 1 IMi'i) 

R. I,. Elman (Brunswick Carp., Costa CA 

9 lh2b) 

(Sponsor: R. P. Cndlr) 

Fine pnrtlrles (froo about 50 to 0.05 micro- 
meter* In diameter) in the eruption plume* from 
two active volcanoes, Ftiego and Sant lagul to, 
were sampled In slut te dele ml nr mas* rwu-cn- 
trulion and site d 1*1 r Unit Ion. with elrn.-ni.il 
composition ami oar pin* logy determined from post 
.sampling analysis. The sample* were taken 
during an aircraft sampling oi«slon In Guatemala 
in February and March I97R. Tho particles were 


emf'lii) vltli ,t 1 P-sfage cascade fn^artot which 
has pie r.00 1 e* trie crystal* ,i» active colUrlMt 
that ptoilnrp, In real tine, electronic signals 
corresponding t« the mass of tutorial deposited 
00 lie- eel lector surfaee. the v.mpte* were 
phot • 'graphed and nnnlvgeil In a scanning election 
mb r«»«t, equipped with an energy dispersive-^ 
x-r.tv analysis, nt ta« hmeni to determine elemental 
compos 1 1 Ion . thr particle* appeared tr be 
divided among two morphological group*. lsr|« 
agglomerates in tlu* large wtre tin t loan and 
less agglomerated part teles In the small site 
Ifflilli'its, Riving tho overall sire distribution 
a hi modal cliarar ter . Mass di-nsltles nf the 
particles were estimated Item their geometric 
sires as seen In the SIM. The large .ltelimrritri 
were c-it Kuril to have mass ■li-i<« 1 1 1 es less than 
rg/en* while the « I rig lets have m.l** density of 
about .V.fc t>(. Supper 1 I ng data I tom nnothei sen- 
*.-r on hoard the nlretall Indicate that rnut of 
the p.irtlele* are liquid 1-rtt1.1l , prcsumaMy 
H (’1 *ud Ild- 1 *4 . Energy Jlsp,-rslve x-rav intlyals *. 
ieve.ils slgnlttcam amounts of Al‘, SI, S, fa and 
('«• with lesser amounts oi Cl, K and It tn the 
large ar< lomerites and Al, St, Fe and Ti tn ths 
subfile run groups. 
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VOLCANIC E (PLOSION PLUMES: 01 NS I TV, TfMPfRA- 
TURt AND PARTICLE CONTENT ESTIMATES FROM PLUME 
MOTIONS 

L. Wt 1 4oi (f nv(ranrx;n(4l Selene*' (k;pt. , Only. 

of Lancaster, Lancaster l A I drR, U.f.) 

S. Self (Popt, of firth St I oncer. , flirtncuth 

Coilegr, Hanover, N. H. 037SS) 

Analyses are reported of 7 eruption clouds 
fron forgo, fllnrjd at one frarr* per second on 
frbruAry ??, I97ft, as part of NCAR's Centra) 
American Volcano Plu^e 5»Tp)lng Experiment. The 
movenirnt of a cloud between frames gives Its 
velocity. The vc loc l ty/t ime curve can be dif- 
ferentiated to find the acceleration and hence 
the tola' force acting on ihe cloud. This force 
consists of 2 parts, air draq (nost Important at 
high speeds early In the nation) and buoyancy 
(most Important late In the motion when the 
cloud is convectlng slowly). Eauatlons of mo- 
tion were derived for 2 cormon cloud shapes, 
sphere' and vertical cylinders, and cloud den- 
sities were calculated by fitting the equations 
to the observed motions. A simplified analysis 
of the heat budget of a cloud permits an esti- 
mate of cloud t*mpe r ature and particle welqht 
fraction to bp m<ida from the density. Clouds 
generally reached temperatures within 10° of 
that of the surrounding air within 10 seconds 
of formation. Dense particle weight fractions 
were less than 2'- by this time, and cloud rise 
velocities less than 5 m/s. In 0 cases it was 
possible to estimate the manimijn siie of dense 
particles supported by the cloud motion: values 
found were 315, 700, HO, and 560 microns. 


v «? I • v 1 ftD ■ 
inset in His in vulcanic imn*niHi clouds 

U, 0, 5nlth 

H. A. /irlln>kl (both at U.S. Geological Survey, 

Denver, Co. '07751 
(Sponsor: C. ». Ilaesrrl 

In I pbruary ami H.»r« b nf |7*s flltrr sampler', 
mnnnlrtl tin an airrraf* »i-rr n’.i-tl In ft»| IrtJ llic 
arr«»5(*l f r.M’l Ion of ll*r crtipllnn clmuli of thtrr 
active Huai rm.il.in volr.mnn (Itirtjn, I’.iraya .pul 
Saul la«p»lto). The filler '.ampler consisted of 
three filters In series. I hr first »a> a I luoro- 
porc filler wilt* ,* nominal pom sir** of (1.5 
microns fur remnv |n«| p. 11 1 1 1 mI.iT r% ( r*** (lie crop- 
I Ion elnnd. Ilig next l wo filler* were Hli.ilnin 
fillers treated with te1ral»,l)| ,rnitil>n« hydrox- 
ide to preferotl lal l> rrlaln aeldlr volatiles 
contained in I he cloud. I|*«* mass of air sar*|*led 
by I he filler.* r.nvj'vl fr«*» . J5 |,» h.(, k»i. lor 
eaeb vole. pm, samples werr collected al JiM-rraslr") 
distances from l he vent. After collect loo of the 
s.injilrs, r.teh filter was rvtrarled with r.o ml of 
water, Sjillts were filtered and .maly/cd fur |9 
n»e1 al S b> Inductively coupled plar-rv*--npl teal 
colssln** speclroractry. 1 luorldc and chloride 
were anal) red l»y specific loo electrode. 

Ihc element* found In significant concentra- 
tions In Ihc 1 ) nor iipore extracts are as follows 
(In decreasing order of cnncenlx at lot*!: Si, la. 

If, M*l. Cl, Ha, I . Mo, /»*, Ph, Co. Sr and Cd. 
These elements represent component* of water- 
soluble sub) invites derived from interaction or 
particulates a»»d acidic gases of the eruption 
cloud, llcornts found In the extracts of the 
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